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Abstract 

In this paper we give a more geometrical formulation of the main 
theorem in [El] on the inverse problem for the second order hyper- 
bolic equation of general form with coefficients independent of the 
time variable. We apply this theorem to the inverse problem for the 
equation of the propagation of light in a moving medium (the Gordon 
equation). Then we study the existence of black and white holes for 
the general hyperbolic and for the Gordon equation and we discuss 
the impact of this phenomenon on the inverse problems. 

1 The refined formulation of the main 

theorem. 

Let Q he a smooth bounded domain in R". Consider a second order 
hyperbolic equation in a cylinder Q x (— oo, +oo): 
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where x = (xi, E Q, Xq G R is the time variable, 

g{x) = (det[(/'''^(x)]"j^,^Q)~^ We assume that g^^{x) = g^^{x) are real-valued 

smooth function in C°°(f2) independent of xq. 
The hyperbolicity of fll.ip means that the quadratic form Yl] k=o 9'^''(^)^J^k 
has one positive and n negative eigenvalues, Vx G fi. We assume in addition 

that 

(1.2) c?°°(x) > 0, X G n, 

i.e. (1, 0, 0) is not a characteristic direction, and that 

n 

(1.3) Yl a'^i^Mk < for V(ei, Cn) 7^ (0, 0), Vx G n, 

j,k=i 

i.e. the quadratic form fll.3p is negative definite. Note that ( 11.31) equivalent 

to the condition that 

(1.4) any direction (0,^) is not characteristic for Vx G il. 

We shall give later another equivalent characterization of the condition 

We shall study the initial-boundary value problem for the equation (11.11) in 

(1.5) w(xo,x) = for Xo << 0, xEfl, 

(1-6) ^lanxR = /, 

where /(xq, x') has a compact support in dQ x R. 
Consider the Dirichlet-to-Neumann (DN) operator 




where i/q = 0, (z/i, z/„) is the unit outward normal to dQ C R" with 
respect to the Euclidean metric, u{xo,x) is the solution of (II. ip . (11.50 . (II. 6p . 

Let Fq C dQ be an open subset of dfl. We shall say that the DN operator 
A is given on Fq x (0, Tq) if we know A/|rox(o,To) for any / with the support 

inrox[0,ro]. 
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Let X be a closed compact subset of dfl x (—00, +00). For each 
distribution / on dQ x (—00, +00), supp f C X, denote by D^{f) the 
support in X (—00, +00) of the solution of the initial-boundary value 
prolem for (11. II) with it = for xq << 0, x G fi, w|anx(-oo,+oo) = /■ 
We define the forward domain of influence of X as the closure of the set 
Usupp fcxD+{f), where the union is taken over all / on dQ x (—00, +00) 

with supports in X. 

We shall give a geometrical description of the set Z)+(X). Let 

[9jk{x)]lk^o = iy%k=o)~\ i-e- [9jk{x)]lk=o is the pseudo-Riemannian 
metric tensor. We shall say that 7(5) = {xq{s) , xi{s) , ...,x„(s)), s G [0,so] is 
a forward time-like path (ray) if 7(5) is continuous and piece-wise smooth 

and each smooth segment 7^*^(5) of 7(3) satisfies 

(1.8) ^^^^^^^^^~ds~ds ^i^^^ ^i+i' 

j,k=0 

^ > 0. Then D+(X) is the closure in x (—00, +00) of the union of all 
forward time-like paths (rays) starting on X. Analogously D^(X) is called 
a domain of dependence of X if D-{X) is the union of all backward 
time-like rays starting at X. We say that piecewise smooth 7(5), s > 0, is a 
backward time-like ray if (11.81) holds for each smooth segment, and ^ < 
on 7^*^(5), Vs G [sj, Sj+i]. Note that the vertical ray Xq = s, x = x^'^\ s > 0, 

is time-like iff gooix^^^) > 0. 

Proposition 1.1. Condition U.3\) holds iff goo{x) > 0, i.e. the ray Xq = 

s, X E Q, s > 0, is time-like. 



Proof: We have 



g ^{x) 



Note that 



j,k=0 



j,k=i 
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Since the quadratic form Yl] k=o d'^'^^J^k has the signature (+1, —1, — 1) 

the quadratic form 



j,k=i 



is negative definite. Therefore Yl]j k=i 9''^i^)^j^k either has only negative 
eigenvalues, or it has one zero eigenvalue, or it has one positive eigenvalue. 

In the first case sgn(det[5'-'^]"j,^;^) = (—1)", in the second case 
det[g^'^{x)]'^i^^^ = 0, and in the third case sgn{det[g^'']'^f^^j^) = (— l)""^. 
Since sgn g{x) = (—1)" we get that goo{x) > iff condition (11.31) holds. □ 

Let r+ be the smallest number such that 
D^(Tq X {xq = 0}) D X {xq = T+}. Analoguously let T_ be the smallest 
number such that D_(Tq x {xq = T_}) D f2 x {xq = 0}. 

We shall require that 

(1.9) To>T+ + T_. 

Consider the following change of variables in f2 x (— oo, +oo): 

(1.10) yo = xo + a{x), 

y = ^{x), 

where y = if{x) is a diffeomorphism of Q onto Q, a{x) G C°°{Q), 

(1.11) ip{x) = / on fo C dQ, 

a{x) =0 on Tq. 

Note that if u{yQ,y) = u{xo,x), where {yo,y) and {xo,x) are related by 
(ll.lOp . (11.111) . then u{yo,y) satisfies an equation of the form 

(1-12) E , ' #- ( VFWWrg^^y)^^^] = 

in n X (— oo, +oo) with initial and boundary conditions 
(1.13) u(yo,y) = for yo«0, yeQ, 

(1-14) «lanxR = /(?/o,l/'), 
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where /(xq, x') = f{yo, y') on Tq x (-oo, +oo) . 
Let A be the DN operator on Tq x (—00, +00) corresponding to fll.l2p . 

dni, (nn. 

The following theorem holds (c.f. [El], Theorem 2.3) 

Theorem 1.1. Let ([TTP, ^LM), (ESi and [riB) . [Tm) . (TJ^ be two hy- 
perbolic initial-boundary value problems in Q x H and Cl x H., respectively, 
and dQndCl D Tq. Suppose conditions U.2\) and U.3\) are satisfied for U.l\} 
and M.12\] . Suppose supp f C Tq x [0,To] in p7^) and (TJ^. If the DN 
operators A and A are equal on Tq x (0, Tq) and ifT^ > T^ + T_ (c.f. 
then there exists a change of variables U.1(J\) . U.ll\) such that 

(1-15) r (y)] -.=0 = J{xW{x)]1,^qJ^{x), 

where yo = xq -\- a{x), y = f{x), 



;i.i6) 



J{x) 



1 a^{x) 



Vx 



is the Jacobi matrix of 1^1. 10\) . 

Note that condition (11. 9p is required for (11. ip only. The condition (11.150 is 

equivalent to 

[9'\y)]-' = {J^{x)r'[g^\x)]-'j-\x), 



or 



;i.l7) [9jki^)]lk=o = J^ix)[g,,iy)]l,^QJix). 

The equality (I1.17P can be rewritten as the equality of differential forms 



:i.i8) 



^ gjk{,x)dxjdxk = ^ gjkiy)dyjdyk. 

j,k=0 j,k=l 



The proof of Theorem II. II is the same as the proof of Theorem 2.3 in [El] 

(see also [E2]). 
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2 The inverse problem for the Gordon 

equation. 

Consider the equation of the propagation of hght in a moving medium. Let 
w{x) = {wi{x) , W2{x) , w^^x)) be the velocity of flow and let 
[v^^\ v^^\ v^'^\ v^^^) be the corresponding four- velocity vector of flow: 

where c is the speed of hght in the vacuum (c.f. [G], [NVV], [LP]). The 
equation for the propagation of light was found in 1923 by Gordon (c.f. 
[G]) and it has the form (11. ip when g^^{x), < j, k < n, n = 3, are 



(2.1) g^\x) = 7f^ + {n^{x) - iy{x)v^{x), 

Vljk\ = [V'^]"^ is the Lorentz metric tensor: rj^'^ = when 
j 7^ k, rp^ = 1, r]^^ = —1, 1 < i < 3, n{x) = e{x)\i{x) is the refraction 
index. Note that [5'jfc(x)] = [fi'-''^]""^ has the form (c.f. [LP]): 

(2.2) 9ik{x) = r]jk + (n"^(x) - l)vjVk, <j,k <n, 

where vq = v^, Vj = —v^, I < j < n. 
The metric (12. 2p is called the Gordon metric and the corresponding 
equation of the form (II. ID is called the Gordon equation. We shall consider 
the Gordon equation in ^2 x (— oo, +oo), where fl has the form 
r2 = ^7o \ U^^ilj, Qq is diffeomorphic to a ball, Qj are smooth domains, 
Qj n fifc = when j ^ k and U^iQj C Q. Domains Qj are called obstacles. 
We shall consider the following initial-bundary value problem for the 

Gordon equation: 

(2.3) u{xo,x) = for xq << 0, x & Q, 

(2.4) u{xo,x)\anjy.-R, = 0, l<j< m, 

(2.5) u{xo,x)\dnox-R. = f{xo,x), 
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i.e. Tq = dflo where Fq is the same as in Theorem ll.il One can consider 
also the case when on some or all obstacles Qj the zero Dirichlet boundary 
condition is replaced by the zero Neumann boundary condition. 
Note that the condition (11.21) always holds for the Gordon equation 

(2.6) (7°° = l + (n2-l)(i;0)2>0. 

The condition (11.31) holds iff 



(2.7) |w(x)p< 



c2 



To prove (12. 7p note that the principal symbol of the Gordon equation is 
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0=0 



Therefore the quadratic form 
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(2.8) -\^\' + in'-l)[j2v'^, 



is negative definite iff (n^ — l)|f P < 1. Since 

v^=^l-\^y\ ^.■= (^1- H!j-^i^^ l<j<n,weget m- 

We shall study the inverse problem for the Gordon equation. The case of 
slowly moving medium (c.f. [LP]) was considered in [El]. The 
Dirichlet-to-Neumann operator corresponding to the Gordon equation has 

the form (c.f. (fLTj) ): 

n " 

A/ = Yl [^'' + ("'(^) - ^^''''^d^.^'^ I - 1 + - i)(E^''^-)'i"^ 

j,k=0 ^ r=l SQoxR 

where u{x) = {ui, Un) is the unit outward normal vector to 

dQo, 1^0 = 0, wlaf^oxR = /• 
We shall impose some restriction on the flow w{x) = {wi{x), Wn{x)). 

Let X = x{s), s G [0, Sq] be a trajectory of the flow, i.e. 
^ = w{x{s)), w{x{s)) ^ 0, s E [0,so]. Consider the trajectories that either 
start and end on OQq or are closed curves in Q. Our assumption is that 

(2.9) The union of all such trajectories is a dense set in Q. 
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Theorem 2.1. Let [5'jfc(a^)]"fc=o ^'^^ [9jk{y)]],k=o Gordon metrics in 

domains VL and Vl, respectively. Consider two initial-boundary value problems 
M.5\) . \1.6\) and U.l^) . U.l^) . Ijl-H^ corresponding to the metrics 
[5'jA:(a^)]",fc=o ^'^'^ [9jk{y)]],k=0' respectively. Assume that that the condition 
^2. 7\ ) holds for both metrics and assume that the refraction coefficients n^{x) 
and n'^{y) are constant. Assume also that the flow w{x) satisfies the condi- 
tion 112. y\) . If the DN operators A and A corresponding to [9'jfc(a^)]",fc=o ^'^'^ 
[9jk{y)]]^k=o ^'^^ equal on dflo x (0,To), where Tq satisfies U.9\) for equation 
U.l\) . then we have 

V? = fi^, Q = Cl, 
and the flows w{x) and w{x) are equal. 

Proof: Analogously to the Remark 2.2 in [El] we can find the symbol of 
the DN operator in the "elliptic" region (c.f. [El]) and retrieve n^{x)\Q^^- 
Since A = A on dflo x (0,To) we get that 'n^{x)\QQ^ = 'n^{x)\QQ^. Therefore 
rt^ = in f2 since we assume that and do not depend on x. Applying 
Theorem 11.11 to the case of the Gordon equation we get that there exists a 
change of variables (11.101) such that a{x) = on dQo, (p{x) = I on dQo and 

ffrrk . ffTTSj) hold where 

(2.10) g^\x) = r]^'' + {n^ - iy{x)v\x), 0<j,k<n, 

~g'\y) = v'" + in' - i)v'{y)v'{y), 0<j,k<n. 

We have from (11.101) 

n 

(2.11) dyo = dxQ + ^a^^{x)dxk, 

k=l 

n 

dyj = ^ ipj^^{x)dxk, 1 <k <n. 

k=l 

Substituting (12. lip into (ll.lSp and taking into account that dxo,dxi, ...,dxn 

are arbitrary we get 

9ooix) = miy)- 

Therefore (c.f. (Q) 

(2.12) 1 + (n-2 - l)vl{x) = 1 + (n"2 _ i)yl^y)_ 
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Note that fll.lSp is equivalent to 



p,r=0 



where ifQ = XQ + a{x). In particular, 



n 



goo^y) = gOO(^x) + 2 ^ c?f°(x)a,^ (x) + ^ [X )a,^a,^[x ) . 

p=l p,r=l 

In the case of Gordon metrics we have 



I + in' - l){d%y)f 

n 

l + {n^- + 2 - l)v\x)vP{x)a,^{x) 

p=i 

n n 

- ^((^x.pi^))'^ +Yi^^- l)v^{x)v''{x)a^^{x)a^^{x). 

p=l p,r=l 



Denote 



n 1 " 

|Va(x)|^ = ^a^^(x), ao{x) = -'Ywp{x)a^^{x). 



c 

p=i p=i 



Since n = n and since fl2.12p implies that f ° = v'^ we get 
ao + 2ao(x)- =0. 

It follows from (12.161) that either 



oo = -1 



|Va| 



K)2(n2-1)' 

or 



Let X = x{s) be a trajectory of the flow, i.e. 

dxk 



Wk{x{s)), 1 < k < n, 0<s<l. 



ds 

Then 

da(x(s)) .^dxk 



^a^,(x(s))— ^ = cao{x{s)), 



ds ^ ds 

k=l 

i.e. cao(x(s)) is the derivative of a{x{s)) along the trajectory of the flow. 
Suppose X = x{s), < s < 1, is a tragectory that starts and ends on Silo- 
Therefore a{x{0)) = a{x{l)) = since a{x) = on OQq. Then ao{x{s)) can 
not satisfy (12.171) since '^"■^^^^^^ = cao{x{s)) < on [0, 1]. Therefore 
^"•^^^^^^ = caQ(x{s)) satisfies the equation fl2.18l) . i.e. '^"■^^^^^^ > on [0, 1]. 
Since a{x{0)) = a(x(l)) = we must have a{x{s)) = and Va(x(s)) = on 
[0, 1]. In the case when x = x{s), < s < 1, is a closed trajectory, i.e. 
x(0) = x{l), a{x{s)) can not again satisfy fl2.17p since '^"■^^^^^^ < on [0, 1] 
and a(x(s)) satisfies fl2.18p only when Va{x{s)) = and a(x(s)) = const on 

[0,1]. 

Since the condition (12. 9p holds and since a\gQ = we get that a = in Q. 
In fact if a(xW) ^ then a(xW) ^ in some neighborhood Uq of x' 
Since Va(x) = in f/g we get that a{x) = a(x'^°^) in Uq, i.e. the set 
a{x) = a{x^^^) is simultaneously open and close. Since a = on O^Iq we get 

that a{x^^^) = 0, and this is a contradiction. 
Applying (12.131) for A; = 0, I < j < n, and for 1 < j, k < n and taking into 

account that a = we get 



n 



(2.19) in' - l)v\y)v^{y) = {n' - l)J2v\x)v^{x)ip,.,^{x), 

p=i 

n 

(2.20) - V'' + in' - l)v'{y)v'{y) = -Y. 

n 

+ ^ (n^ - l)vP{x)v''{x)(fj^^ipkxr, l<j,k<n. 

p,r=l 

Using again that n = n and = we get from (12.191) 

n 

(2.21) v^{y) = J2v''{x)^j,,{x), l<3<n, d\y) = v\x). 

p=i 
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Substituting fl^:^ into fl^T^ we obtain 



(2.22) rf^ = Y,^,xM^k-M^ l<j,k<n. 

p=i 

In particular, we have 

n 

(2.23) ^%^^) = 1' "Pjldno = Xj, I < j < n. 
p=i 

The Cauchy problem fl2.23p has a unique solution in 
Q : (fij = Xj, 1 < j < n. Therefore y = ip{x) = x in f2. This implies that 
fl = Q. Also (pj{x) = Xj implies (c.f. (12.211) ) that v^{x) = v^{x), I < j < n. 
Therefore Wj{x) = Wj{x), 1 < j < n. □ 

3 Optical black holes. 

In this section we explore the situation when the condition (11.31) is not 

satisfied. 

Denote by S the surface in Q given by the equation A(a;) =0 where 

(3.1) Aix) = dety\x)]l,^,. 

We assume that is a smooth and closed surface. Denote by flext the 
exterior of 5" and by flint the interior of S. We assume that A > in 
flext n f2, A < in flint near S. It follows from Proposition II. II that the 
equation of S can be written in the form goo{x) = and goo{x) > in 
flext n fl, goo < in flint near gooi^) = 0. We shall write often the equation 
of S in the form S{x) = and we assume that is an outward normal 

to the surface S{x) = 0. 
In the case of Gordon metric the equation of S has the form 

n 

{n\x)-l)Y,iv'i^)?-^ = 0, 



or, equivalently. 



j=i ^ > 
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and the domain Qint is given by the inequahty 

\ w'^[x) > 

■' n^ix) 

j=i 

near S{x) = 0. 

Suppose that S{x) = is a characteristic surface of the equation (11 .ip . i.e. 

n 

(3.3) J29''i^)S.,{^)Sx,{x) = when S{x)=0. 
j,k=i 

Let Xj = Xj{s), = ^j{s), < j < n, s > be a null-bicharacteristic of 

dLH), i.e. 

(3.4) ^ = 2j2g'\x{s)C,{s), x,(0)=y,, 0<j<n, 

k=0 



(3.5) ^ = -J29i'M'^)U-^^Us) m=Vv 0<p<n, 

j,k=0 

where x{s) = (xi(s), ...,x„(s)), r^o = 0, t] ^- The null-bicharacteristic 

means that 

n 

(3.6) 5^ 9'\x{sm,{s)Us) = 0. 

Note that if 

n 

(3.7) 5^ 9'\y)VjVk = 

then ([Ml) holds for all s G R. 
The bicharacteristic (null-bicharacteristic) is a curve in 
To(R"+^) = R"+^ X (R"+i \ {0}) and its projection on R"+i is called a 
geodesic (null-geodesic). We shall show now that the null-geodesic satisfies 

the equation: 

(3.8) ± f ^ 0. 

j,k=0 
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Denote by G the matrix [g^^{x)]^^f^^Q. Then G ^ = [gjk{x)]^^i^^Q. It follows 
from ([S3D that f = 2G'f(s), where x{s) = {xo{s),x{s)), (= {^o,0- Then 

^"'= iG-ig and (EJ]) implies 

that is equavalent to (13. Sp . 
Fix an arbitrary point y on ^(x) = 0. Denote by K^{y) the following 

half-cone in R""*"^: 

(3.10) K^{y) = {e = (^0,0 e R"+^ : {G{y)tO > 0, > 0}. 

Let K^{y) be the dual half-cone in R"~*"^: 

(3.11) K+{y) = {x= (±0, x) e R"+^ : {G~^x, x) > 0, io > 0}. 

Note that 

(3.12) (i',f)>0 for any x e K+{y) and any ^eK+{y). 

Consider the nuU-bicharacteristic {x{s),^{s)) (c.f. (13. 4p . (13. 5p . (13.61) ) with 

the following initial conditions: 

(3.13) xo{0)=yo, x{0)=y, Co{0)=0, m = SM, 

where S{y) = 0. We have (c.f. ^^): 

(3.14) '^ = 2±g'^^iy)S.M 
Proposition 3.1. Let {0,b{y)) be a characteristic direction, i.e. 

n 

Y,9'\y)bAy)h{y) = o. 

Suppose 

n 

Then (0, b(y)) G K^{y) and the half-cone K+(y) is contained in the half-space 
{(ao,ai, ...,a„) : ZlJ=i > 0}- 
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Proof: Let e > be small. Then {e,bi{y),b2{y), ■■■,bn{y)) G K^{y) since 
g^%y)6' + 2eY:,=^9'\yMy) > (c.f. dSZUD). For any (xo^x) E K+{y) we 

have (c.f. (1312]) ): 

n 

xqs + ^Xjh,j{y) > 0, 
j=i 

i.e. K^{y) is contained in the half-space eaQ + YTj=i '^j^jiv) — 0- Taking 

the limit when e — > we prove Proposition 13.11 
It follows from the hyperbolicity of (11. II) and from (13. 3p that the right hand 
side of (13.141) is not zero. In fact, the hyperbolicity implies that the 
equation Y^j k=o 9''''(y)^j^k = has two distinct real roots ^o'\C),^(^\0 ^r 
any ^ ^ 0. Taking ^ = S^{y) we get ^oo(^)^2 ^ 3 9'Ky)^oS-,^ iv) = 0, 

i.e. = 0, = -2{9''{y))-'EU9''iy)s.M ^0- ° 

We assumed that S^ix) is the outward normal to the S{x) = 0. 
Since S{x) satisfies (13. 3p when S{x) = we have that either 

n 

(3.15) J29''iy)S^M>^^ such that S{y) = 0. 

or 

n 

(3.16) Y.9''iy)S.M<^^ S{y) = 0. 

i=i 

If (13151) holds then by Proposition O (0, S^{y)) E K+{y), My, S{y) = 0, 
and K+{y) is contained in the half-space (0, Sx{y)) ■ {xq, x) > 0. 
All directions of Kj^iy) except (1, 0, 0) are pointed inside VLext x R- 
Therefore the forward domain of infiuence of the surface {>S'(a;) = 0} x R 
does not intersect f2j„t x R. In such case the surface {S{x) = 0} x R is 
called the boundary of a white hole. Note that the forward domain of 

infiuence of Vt-int x R does intersect ^ext x R- 
Consider now the case when (I3.16P holds. Then (c.f. Proposition 13.11) the 
dual half-cone K+{y) is contained in the half-space x ■ {0, —Sx{y)) > and 
all X E K^{y) (except (1,0, ...,0)) are pointed inside VLmt x R. Therefore 
the domain of infiuence of VLint x R is contained in VLint x R. 
The surface {S{x) = 0} x R is called the boundary of a black hole in this 

case. 

We proved the following theorem: 
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Theorem 3.1. Let S{x) = det^'']],^^^ = is a closed and smooth charac- 
teristic surface, goo{x) > in Qe^t H Q, goo{x) < in Qmt near S{x) = 0. 
Let {0,Sx{x)) be the outward normal to the surface {S{x) = 0} x R. Then 
{S{x) = 0} X R forms a white hole if Ii3.15\) holds, and a black hole if Ii3.16\) 
holds. 

The apphcation of Theorem 13.11 to the Gordon equation yields the following 

result (c.f. [V]): 

Theorem 3.2. Let S{x) = be the surface X]j=i'"^|(^) = Ti^x) 
{S{x) = 0} X H. be a characteristic surface for the Gordon equation. Then 
{S{x) = 0} X R forms a white hole if w{x) = {wi{x), ...,Wn{x)) is pointed 
inside Qext when S{x) = and {S{x) = 0} x R forms a black hole if w{x) is 
pointed inside Qint when S{x) = 0. 

Proof: Since {S{x) = 0} x R is a characteristic surface for the Gordon 

equation we have (c.f. (12. 8p ): 



X) 



Therefore 



(3.17) (c^ - \wnSx\' = in' - l)iwix) ■ 5,.(x))^ 



Since = when S{x) = we get IS^^plwp = (w ■ Sx^. This last 
equality holds iff Sx{x) = a{x)w{x) where S{x) = 0, a{x) ^ 0. Since Sx{x^ 
is an outward normal we have that a[x) > if w[x) pointed outwardly, and 
a{x) < if w{x) is pointed inwardly. In the case of Gordon equation 

^oi(x) = f'^-^l Wj{x) (c.f. ^Mi). We have 

'YTj=\9^\^)^xAx) = a{x)\w\''. Therefore (I3.15P holds when w{x) is 



pointed into Qext and (13.161) holds when w{x) is pointed into flint- In the 
first case we we have a white hole and in the second case we have a black 
hole. □ 
We shall study now the impact of the existence of the black or white hole 

on the uniqueness of the inverse problem. 
Consider the case of a white hole. Then the domain of the dependence (i.e. 
the backward domain of influence) of any point {yo, y) G VLmt x (— oo, +oo) 
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is contained in flint x (—00, +00). If u{xo,x) is the solution of the 
initial-boundary value problem fll.ll) . (11.51) . (11.61) we get that -u = in 
^int X (— cxo, +cxd) by the uniqueness of the Cauchy problem. Therefore the 
change of the coefficients of (II. ip in Qint x (— cxd, -|-cxd) does not change the 
solution of the initial-boundary value problem (11.11) . (11.51) . (11.61) in 
{Q n flext) X {—00, +00) and it does not change the Cauchy data on 
OQq X (—00, +00). Therefore the boundary data on dflo x (— cxo, +00) do 
not able to determine the coefficients of (II. ip . (modulo (ll.lOp . (11.111) ) in 

flint, i-e. we have a nonuniqueness. 
Consider now the case of a black hole. Then the domain of influence of each 
point (2/0, y) G flint X R is contained in flint x R and the domain of 
dependence of each point {yo,y) G f^ext x R is contained in flg^t x R. 
Consider the initial-boundary value problem (II. 5p . (11.61) for two hyperbolic 
equations L'^^^mW = in f2 x {—00, +00) whose coefficients differ in ^int- 
Since domain of dependence of any (j/o, y) G flext x R is contained in 
ilext X R the solutions u^'^\i = 1,2, of the initial-boundary value problem 
fll.51) . (11.61) are equal in fl^xt x {—00, +00). Therefore the A'^^) and A*^^) are 
equal on dflo x {—00, +00), i.e. the boundary measurements are equal 
despite that the coefficients of L^^^ and L^'^^ differ in ^lint- Therefore we 
again have a nonuniqueness of the solution of the inverse problem. 
Change the formulation of the initial-boundary value problem allowing 
nonzero initial condition in fl^nt x R or allowing the nonzero right hand of 
(II. ip with support in x R. Then u{xq, x) will be nonzero in 
f^int X (—00, +00) but this will not effect the solution in flg^t x R. 

4 Black holes inside the ergosphere. 

Let S{x) = 0, flint, ^ext be the same as in §3. Borrowing the terminology 
from the general relativity we shall call S{x) = the ergosphere and 
fie = ^int n {A(a:) < 0} n f2 is ergoregion (c.f. [V]). We assume that 

A(x) <Omfle\S. 

If 5" = {x : S{x) = 0} is a characteristic surface then {flint fl f2) x R is 
either a black or white hole (c.f. Theorem 13. ip . Now consider the case when 

S is not characteristic at any point y E S, i.e. 

n 

j,k=i 
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where u{y) is the outward normal to S. Note that the quadratic form 
YTj k=i9^^^-'^)^3^k has the signature (1, —1, —1) when x G Vt^. The 
question is whether there exists a black or white hole inside VL^ x R. 
Consider the case n = 2, i.e. the case of two space dimensions In this case 
there are two families of characteristic curves S"^ = const: 



(4.2) Yl a^'i^K^. 



± 

j,k=l 



in a neighborhood of any point of Qe- Equation (14.21) can be factored and 

we get in the region = {x G fie : 9^'^{x) ^ 0}: 



Note that A(x) = g^^{x)g'^'^{x) - {g^'^{x)f < in Vt^- 
We shall derive an equation for the characteristic cirves S^{x) = const that 
holds in a neighborhood of any point in fie- Let Ui be the set in fie where 
either g'^^(x) 7^ or g^'^{x) > 0. Denote by f^^{x) a nonzero vector field 
fu^{x) = {g^'^{x) + ^y—A{x),g'^'^{x)), x e Ui. Let U2 be the set where 
either g^^{x) 7^ or g^'^{x) < and denote by f^^{x) the nonzero vector 
field /+ (a;) = {g^\x) , g^^ (x) - V^A(^), x G U2. Let 
f/f = fie \ Ui, i = l,2. Note that n UI = {x : c/" = g^^ = g^^ = 0} = 
since the rank of [g-''']'jk=i is at least 1. Therefore U1UU2 = fie- Let 
Vi, i = 1, 2, be a small neighborhood of t/f , i = 1,2, such that VinV2 = 0- 
Denote Ui = Ui\ Vi, i = 1,2. Then f/^ U f/a = fie- Note that 



f ^ g^^-V^ \f \Un ■ T~T rM~T 

g^^{x) + y/-A g^\x) 



since ^"(?22^= - ^f^W^ + v^)- 
Extend X{x) as nonzero function from Ui fl f/a to such that A(x) is 
continuous in U2 and smooth in U2 \ S, and define 

(4.4) /+(x) = /+(x) in [/i, /+(x)=A-i(x)/+(x) in U2. 

Then /"*" is a continuous nonzero vector field in fie that is smooth in fie \ S. 
Analogously let /^^(x) = {g^^ - g''),fuS^) = ig'\ g'' + V^), 
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where is the set where either g^"^ 7^ or g^"^ < and U4 is the set where 
either g^^ ^ Q or g^^ > 0. Then fj}{x) = Xi{x)f^{x) in U3 f] U4 and Ai ^ 
in U3 n U4. Here U3 C U3, U4C are similar to Ui, i = 1,2, in fj44|) . 
Extending Ai (x) from f/3 fl f/4 to [/4 we get a vector field /~ {x) in f2e that 

is continuous in fig and smooth in \ S. 
We have that fl4.2p is equivalent to 

(4.5) ft{x)S^^{x) + fHx)Si{x)=0, 

where /^(x) = {f^{x), f2{x)). Note that /^(a:) 7^ (0,0) for any x G fig and 

/+(x)^r(x),Vxefie\>S. 

Since the rank of [g^'^ly)]^ j^^i is 1 on there exists a smooth and nonzero 

Ky) = {bi{y),h2{y)) such that 

2 

(4.6) Y.9'\yMy) = ^, J = 1,2, yes. 

k=l 

Since 

2 

(4.7) 5^<?^^y)6,(2/)6,.(|/)=0 

j,k=i 

we get from (14. II) that b{y) is not colinear with ly^y), G S. It follows 

from (14. 5p with S^{y) replaced by b{y) that 

(4.8) ft{yMy) + fHyMy) = o, ^y e s. 

Note that /^(y) = f^{y) when y e S. Therefore the condition (14. ip implies 
that f^{y) is not tangential to S for any y e S. Therefore changing f^{x) 
to —f^{y) if needed we will assume that f^{y) is pointed inside Qe for any 

ye 5. 

Consider differential equations 

(4.9) ^^ = /±(a;±(a)), x^{0) = y, y e S, a>0. 

da 

Let S'^(x) = const be a characteristic curve. Then fl4.5p . fl4.9p imply that 
£S^{x^{a)) = 0, i.e. S^{x^{a)) = S^{y),a> 0. Therefore x = x^{a) are 
the parametric equations of the characteristic curves S^{x) = S^{y). It was 
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shown in §3 that X]j=i 7^ if b{x) = {pi, 1)2) satisfies fl4.7p . 

Changing h{y) to —h{y) if needed we assume that 

2 

(4.10) Y.3''iyMy)>^^ VyG5. 

We shall denote by 5"^ x R the characteristic surfaces in fig x R where 
S^{x) = const and xq G R. Let {yo, y^) be any point of x R. Consider 
the nuU-bicharacteristic (13 ■4p . (I3.5P with initial conditions 
4(0) = 2/o,x±(0) = ?/±,^±(0) = St{y^).i^{Q) = 0. Note that 
(s), s > is the corresponding null-geodesic. We shall 
show that S^{x^{s)) = 5*^(2/), i.e. this null-geodesic remains on the 

characteristic surface x R. We have 

2 

since T,lk=o9''vfVk = 0, = 0, = S^{y^) (c.f. dSZD). Note 

that ^q{s) = C^{0) = since [g^''{x)]'^j^i^^Q is independent of xq (c.f. f l33|l ). 

Since S^{x) is a solution of the eiconal equation 
Yl'j k=i9'^^i^)'^xji^)Stki^) = in an open neighborhood of in fig, we 

have that (c.f. [CH]) 

CHs) = S^{x^{s)). 

We have (^^, = 2G{0,^^{s)). Therefore 

x: ^s^M^is)) = f: ^^fis) = 2 x: ^^^^i-^i^))?? (^)e.(^) = o, 

j=l j=l j,k=0 

i.e. f 5±(x±(s)) = 0, Vs. 
Therefore the projection of the null-bicharacteristic 
Xq = Xq{s),x = x^(s), x^(0) = y^, Co = 0, ^ = ^^(s) ou the {xi, X2)-plane 

belongs to the curve S'^(x) = 5^(1/^). 
Fix arbitrary y E S. Denote by Il~^{y) the intersection of the half-plane 

(4.12) CMy) + My) > 
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with fig- Analogously let 11 (y) be the intersection of the half-plane 

(4.13) CMy) + My)<o 

with Qe- Let x = a;^(cr) be the solution of (jH]), £^(0) = y,a>0. Note 
that f~^{y) = f~{y) G n^(y) n n {y). However for cr > and small, either 
x~^{a) belongs to 11+ (y) and x~{a) to Il~{y), or vice versa. For the 
definiteness let x^{a) G Il'^{y) for < a < e and x~{a) belongs to n~(y) 
for < a < e. Condition (14.11) and the continuity in y imply 
x+(cr) e n+(y), x'{a) G Il~{y), < a < e, for all y G S". Each x=^(cr) is a 
parametric equation of the characteristic curve 5'+(a;) = S'^^y), i.e. 

(4.14) 5^(£^(a)) = 5^(2/), Va>0. 

Note that the curve S'^{x) = S'^{y) is contained in 11+ (i/) when x ^ y and 

close to y. 

Consider now two null-bicharacteristics 

(4.15) xo = x^is), x = x+(s), eo = 0, C = tis),s>0, 

with initial conditions Xq{0) = 0,a;+(0) = y,^^{0) = b and 

(4.15') xo = Xq{s), x = x~{s), ^0 = 0, ^ = C{s),s>0, 

with initial conditions Xg (0) = 0,x^(0) = y,^~(0) = —6. 

Note that 

(4.16) ^^ = 2j2g'%x^smi{s)>0 



j=i 

since fl4.10p holds. Also we have that 



(4.17) ^^<o, s>0, 

as 

since ^~(0) = —b. Therefore Xq = Xq{s), x = x~^{s) is a forward 
nuU-geodesics since xq is increasing when s is increasing and 
X = Xq{s), X = x~{s) is a backward nuU-geodesics. 

It follows from ([33]) and (gl]) that 

=2 5^^?^'^(y)6.(i/) = 0, J = 1,2. 



ds 

k=l 
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However 

^^ = 2G{x^{s)){0,e{s))j^0, s>0, 
as 

since A{x) 7^ in f2e \ S. 
Since {0,b{y)) G K^{y) (c.f. fl3.10p ) the dual half-cone K^{y) is contained 
in n {y). Therefore the projection x = x'^{s)^ s > 0, of the 

nuU-bicharacteristic fl4.15l) satisfies 

(4.18) S+{x+{s)) = S+{y), s > 0. 

Analogously, the projection x = x~{s) of the null-bicharacteristic (4.15') on 

(xi, X2) -plane satisfies 

(4.18') S-{x-{s)) = S-{y), s>0. 

Comparing fl4.14p and (14.181) . (4.18') we get that x = x'^{s),x = x^((t) are 

different parametrization of the same curve, i.e. 

(4.19) x^{s^{a)) = x^{a), a > 0, 

where s±(0) = 0, > for (x > 0. □ 

Remark 4.1 Note that if 

(■5)5^0 = 0,^ = ^(s), satisfies 
( 13. 4p . (13. Sp . then x = xq{—s),x = x(— s),^o = 0,^ = — 1^(— s) also satisfies 
(13.41) . (13.51) . Therefore changing s > to — s in (4.15') and combining (14.151) 

and (4.15') we get a forward null-bicharacteristic 
X = xo{s),x = x{s),^o = 0,^ = (^(s), defined on (— cxd, +00) with initial 
conditions xo(0) = 0,x(0) = y,C.{0) = b. The projection of this 
bicharacteristic on (a;i, a;2)-plane has a singularity (caustic) at x = y. □ 
Let Qe C Qn Qint n {A < 0}. We assume that the boundary of Qe consists 
of = dQint and Si. Assume that A < on 5*1 and fie is diffeomorphic to 

an annulus domain in R^. 
Let y be an arbitrary point on Si. Consider the forward cone of influence 
K^{y). Let K^\y) be the projection of this cone on the (xi, X2)-plane. We 

assume that 

(4.20) N{y) ■x>0, Vy G ^i, Vi; G (y), 

where N{y) is the outward unit normal to 5*1 and x = (xi,i;2) ^ is any 
vector in Kf\y). In particular. Si is not a characteristic curve. Note that 
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condition fl4.20l) is equivalent to the condition that 



(4.21) N{y)-^^^>0, x^{s^)=y,yyeS„ 

where x^{s) is the projection on (xi, X2)-plane of two forward 
nuU-bicharacteristics such that x = x^{s) are parametric equations of two 
characteristics of the form (14. 2p passing through y & Si. Note that 

^^4?^ > when x±(sf ) G ^i. 
We shall say that conditions (4.20') or (4.21') are satisfied if N{y) ■ i: < in 

(g^DD or N{y) ■ < hold in flOT]l . 

Theorem 4.1. Suppose dQe = 'S'o U Si, where A(x) = on S, A(x) < on 
Qe \ S. Suppose ( |y^.i[ ) holds on S. Suppose also that either ( [^.^0[ ) or (4-20' ) 



is satisfied on Si. Then there exists a Jordan curve So{x) = between S and 
Si such that So x H is a characteristic surface, i.e. Sq x H is a boundary of 
either a black or a white hole. 

Proof: It was shown already that condition (14.11) implies the existence of 
two null-bicharacteristics Xq = Xq (s),x = x^(s),^^ = 0, ^ = ^^(s), s > 0, 
such that x^{0) = y & S and x^{s) after reparametrization 
s = s^((t), s^(0) = 0, "^^Jf^ > 0,a > 0, coincide with the solution x = x^((t) 
of the differential equation (14. 9p : x^(s(cr)) = x'^{a),a > 0. Moreover 

^" < 0, i.e. xq decreases when a (or s) increases, and > 0. 

Suppose condition (I4.20p is satisfied. We shall show that there is no 

solution X = x~{a) of 

(4.22) = /-(£(.)), 

da 

x{0) = y E S that reaches Si, i.e. x'{ai) = y^^^ where y^-^^ G Si,ai > 0. 
Suppose such x = x~ (a) exists. When a > ai x = x^ (a) leaves fl^ since 
Si is not characteristic. Note that for the nuU-bicharacteristic whose 
projection is x = x~{s{a)) = x{a), the time variable xq = Xq (s) is 
decreasing when a is increasing. Therefore x = x~{a) leaves Qe when xq is 
decreasing. From other side the condition (14.201) implies that the projection 
of all null bi-characteristics passing through y^^^ G Si leave fig when xq 
increases. This contradiction proves that the limit set of the trajectory 
X = x~{a) is contained inside Qe- By the Poincare-Bendixson theorem (c.f. 
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[H]) there exists a limit cycle, i.e. a closed periodic solution x = z~{a) of 
f l4.22p which has no points of self- intersection. Let So{x) = be the 
equation of this orbit. Then 5*0 x R be a characteristic surface where 
So = {x : 5*0 (x) = 0}. Other solutios of (14.221) are spiraled around 5*0 when 

xo infty. 

Now we shall assume that the condition (4.20') (or (4.21')) is satisfied. 
These conditions mean that the projection on (xi, a;2)-plane of any 
null-bicharacteristic passing through 5*1 enters Qe when xq increases. From 

other side, consider the solution of 

(4.23) ^^ = r{S:H^)), cr>0, 

that starts on S : x(0) = y E S. We claim that this solution can not reach 
5*1. In fact if X = x~^{a) reaches 5*1 it will leave Si when a > ai. We have 

> for the null-bicharacteristics whose projection on the 
(xi, X2)-plane is x = x~^{a) after the change of the parameter s = s~^{a). 
Therefore x = x~^{s) leaves fie when xq increases. This contradiction shows 

that the limit set of x = x~^{a) is contained inside fie- 
Again by the Poincare-Bendixson theorem there exists a closed 
characteristic curve 5*0 (x) = without points of self- intersection. All other 

solutions of (I4.22P spiral around 5*0 when xq +oo. 
Remark 4.2 Without additional assumptions it is impossible to specify 
whether 5*0 x R is a boundary of a white or a black hole. 
For example, assume that the condition (I4.20p is satisfied. Then 5*0 is a 
limit cycle for the equation ^^-j^ = /~(x~(cr)). Since /~(x) ^ /"""(x) for all 
X e we have that /+ (x) , x G 5*0 is pointed either into the interior of 5*0 
for all X e S'o or into the exterior of 5*0. In the first case 5*0 x R forms a 
black hole and in the second case S'o x R forms a white hole. In the second 
case we can apply the Poincare-Bendixson theorem to the equation (14.231) 
in the annulus domain between S and S'o. Then there exists a limit cycle S2 
for the equation '^^^j'^'' = /+(x(cr)) between S and 5*1 and 5*2 x R is the 
boundary of a white hole. Also there exists a white or black hole formed by 
^3 X R where S'3(x) = is a limit cycle for (14.231) between S'o and Si. 
Analoguously if condition (4.20') is satisfied then the limit cycle S'o for the 
equation ( 14.23^ forms a white hole if /~(x) restricted to S'o is pointed in the 
exterior of S'o and there are two additional black or white holes if f~{x)\so 

is pointed into the interior of S'o. 
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Remark 4.3 The black or white holes described in Theorem 14. II are stable 
in the sense that if we perturb slightly [9'"''^(a^)]j,fc=o then assumptions of 
Theorem 14.11 remain valid and therefore the limit cycle solution will still 
exist. This is in contrast with black or white holes obtained by Theorem 
I3.lt when we perturb [g-''']^^k=o then the surface S = {A(x) = 0} may cease 
to be characteristic and the black or white hole at x R disappears. 

□ 

Now we shall formulate the application of Theorem 14.11 to the case of the 

Gordon equation. 

Theorem 4.2. Let dile = S Li Si. Assume that \w{x)\'^ = ^^|^ on S and 
\w{x)\'^ > in ^lf,\S. Suppose w{y) is not colinear with p{y) on S, where 
v{y) is the outward unit normal to S. Suppose that either 

(4.24) - l)^(^;(x) ■ iV) > 1 on 5i, 

or 

(4.24') {n^{x) -iy^{v{x) ■ N) < -I on S^ 

where v{x) = ^1 — ^"'[.2''^ j ^ ^^^5 N{x) is the outward unit normal to Si. 

Then there exists a Jordan curve So{x) = between S and Si such that 
So xH is a characteristic surface, i. e. 5*0 x R is the boundary of either black 
or white hole. 

Proof: Note that A(y) = is equivalent to l^op = ^1^; y G S in the case 
of the Gordon equation and A(x) < implies > ;^^|^ in il^ \ S. 
Since A(y) = there exists a smooth b{y) 7^ 0, y E S, such that 
Yl^k=i9^^iy)^k{y) =0, j = l,2(c.f. fl4.6p ). For the Gordon equation we have 

-hj{y) + {n^ -l){vb)v, = Q, J = 1,2, 

i.e. b{y) is colinear with w{y) on S. Therefore h{y) is not colinear with v{y) 
since w{y) is not colinear, Vy G 5*. Note that condition fl4.10p has the form 
[n? — l)(v ■ h)v^ > 0, i.e. b and w has the same direction. We take b = w. 
Therefore f'^{y) = f^{y) is orthogonal to w{y) and is not tangential to 

S, Wy e S. 
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Let y^^^ be any point on 5*1 . Suppose fl4.24p holds. Consider any 
null-bicharacteristic whose projection is passing through y^^\ The equation 



2 

j,k=i 

has the following form for the Gordon equation 

(4.25) - + - ■ Civ^W = 0- 

Note that = 2{n^ - \)v^{v ■ > implies that {y{x{s)) ■ i[s)) > 0. It 

follows from fH725l) that 



(4.26) V ■ e 



(n2-l)5 

We have 



(4.27) 

dx. 
ds 



"^^^^ - 2^29'% = -2^j+2{n^-l){vOv, = -2^,+2{n^-l)^\^\v„ j = 1,2. 



k=l 



Take the inner product of fl4.27p with A^, where is outward unit normal 

to Si. We get 

^■N = -2^-N + 2|e|(n2 - 1)^(^; ■ A^). 
Since |^ ■ < |^| we get, using (S21, that ^ ■ > on ^i. Note that 



dxo(s) 
ds 



> 0. Therefore any null-bicharacteristic escapes fie when s > Sq (or 

when Xq > Xq^Sq)) where x{so) G 5*1. 
Therefore by Theorem 14.11 we have a characteristic surface 5*0 x R. If 
(n^ — l)2(u • A^) < —1 then all nuU-bicharacteristics are entering when 
xq increases. Therefore again we can apply Theorem 14.11 
As in Remark 4.3 the black and white holes described in Theorem 14. 2 1 are 

stable. 

Now we shall refine Theorem 14.21 and specify when the characteristic 
surface 5*0 x R is the boundary of a white hole and when it is the boundary 

of a black hole. 
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We shall say that the flow = {w, {x),W2{x)) is incoming if for any 
closed simple curve F C fig containing Si there exists at least one point 
y eT such that w{y) ■ v{y) > where iy{y) is a normal to F pointed inside 
F. Analogously, the flow w{x) is called outgoing in Qe if for any F there 

exists ?/ e F such that w{y) ■ iy{y) < 0. 

Theorem 4.3. Consider the Gordon equation in fig where is the same 
as in Theorems ^.i and \4.2\ Suppose w{y) is not colinear with the normal 



to S, ^y G S, and ( [^ . 24^ holds on Si. Assume in addition that the flow 
w{x) is incoming. Then there exists a black hole bounded by Sq x H, where 
Sq is a Jordan curve between S and Si. If (4.24') holds and the flow w{x) is 
outgoing then there exists a white hole with the boundary Sq x R. 

Proof: Suppose (14 .24^ holds. It was proven in Theorem 14.21 that there 
exists a characteristic surface 5*0 x R. Now using that the flow w{x) is 
incoming we shall prove that iSq x R is the boundary of a black hole. 
Take any y e Sq. Let bi{y) be a normal to 5*0. Choose the direction of bi{y) 
such that EUa'^iyMy) = 2(^'(?/) - l)v%y)iviy) ■ bi{y)) > 0, i.e. 
bi{y) ■ w{y) > 0). It follows from the Proposition 13. II that K^{y) (c.f. 
(13.111) ) is contained in the half-space {ao,ai,a2) ■ (0, 62(2/)) > 0. 
Since the flow w{x) is incoming there exists yo G 5*0 such that 
w{yo) ■ iy{yo) > where iy{yo) is pointed inside Sq. Since bi{yo) ■ w{yo) > 
and bi{yo) is colinear with ly^yo) we get that 6i(?/o) is also pointed inside 5*0. 
Since bi{y) is continuous in y we have that bi{y) is pointed inside Sq for all 
y G 5*0. Therefore {xq,x) G K^{y) are pointed inside Sq x R, i.e. 5*0 x R is 

the boundary of a black hole. 
Suppose that (4.24') is satisfled and the flow w{x) is outgoing. Let biiy) be 
the same as above, i.e. hi{y) is a normal to Sq and bi{y) ■ w{y) > 0. Since 
the flow w{x) is outgoing there exists yo G Sq such that ty(?/o) " ^{yo) < 
where z/(|/o) is pointed inside 5*0. Then 6i(|/o) is a normal to 5*0 that is 
pointed outside of Sq. Therefore K^^yo) is pointed outside of Sq x R. The 
same is true for any y E Sq hj the continuity. Therefore 6*0 x R is the 
boundary of a white hole. □ 

Example 4.1 (Acoustic black hole (c.f. [V])). 
Consider a fluid flow in a vortex with the velocity fleld 

A B - 

(4.28) V = {v^.v"^) = -f + —e, 
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where r = \x\, f = , 9 = , A and B are constants. The 

inverse of the metric tensor in this case has the form 



(4.29) = g'^ = g^' = -v^ , I < J < 2, 

pc pc 

^^■' = -(-c% + t;V), l<3,k<2, 
pc 

where c is the sound speed, p is the density. 
It is convenient to use polar coordinates in (14.21) . Assuming c = l,p = 1 we 

have 

or (c.f. (USD) 



(4.31) (-^-ll^W^i '^'^^^ '^''^ 



for r < \/ + B"^. The ergosphere, i.e. the curve where A = has the form 
A? + B"^ = r"^. We shall consider the domain < r < + B'^ or any 
compact domain fig = {ri < r < tq} where tq = + -B^, n < |v4|. The 
condition (14. ip is satisfied when B ^ 0. The condition (I4.20p is satisfied 
when A > ri and the condition (4.20') holds when A < —ri. 
Consider the case ^4 > 0, S > 0. The system (14.91) has the following form in 

polar coordinates: 

dr+ ,0 o d9+ AB 



(4.32) — = A^-r\ — = + + B^ 



ds ds 



r 



2 



d^^ de-{s) _ 1-^ 

^ ds ' ds ^ + Vyl2 + 52-r2' 

r±(0) = ro, ^^±(0) = ^, = tq, s > 0. Note that f = f f + rf ^. We 
wanted to ensure that the vector field in the right hand sides of (I4.32p . 
(14:3311 is not zero. Since ^ - ^JWTW^^ = when A = r, S > we 

divided by — A"^ in (14.33!) and used that 



27 



AR _ + -r"^ = AB_^^J^J^;'2- It is clear that = is the hmit 

cycle for the equation f l4.32p . 
In the notations of Theorem SiKc.f. (I410D . f09|) ) we have 
b{y) = v{y) = ^f + where \y\ = tq = ^W+W. We have, (c.f. (|412|) . 

611): 



as \ as as 









He) 




To J 



Byr^ - + 0(ro - r(s)) > for < s < £. 

dx~ (s) 



Analogously • b{y) < for < s < e, e is small. Note that 

r = r+(s), 9 = 9~^{s) is the projection on (xi, X2)-plane (after a 
reparametrization) of the null-bicharacteristic with initial conditions 

x+(0) = y, \y\ = ro,Xo+(0) = 0,^0^(0) = 0,^(0) = b{y) and ^ > 0, Vs > 0. 
When Xo —>■ +oo, r = r^{s),9 = 0^{s) spirals toward the limit cycle 
r = A. Also r = r^{s),9 = 9^{s) is the projection (after a 
reparametrization) of the null-bicharacteristic with initial conditions 

x-{0) = 0,x-(0) = y,^o = 0,^(0) = -b{y) and < 0. Note that 

r = r~{s),9 = 9~{s) reaches r = ri when xq — oo. Therefore r = A is the 

boundary of a white hole. 
When A< -ri, B > we have, analogously to iK32^ . K3^ : 

(434) '^^-1 '-^ 



(ir (s) ,9 r, d9 is) AB /—- 

(4.35) ' = A^ - r^, — = + -r2, 



(is (is 

(0) = To, 6*^(0) = ill , = To, and we modified (I4.34p since 
M. + _^ _ ^2 ^ Q ^j^gj^ ^ ^ 1^1^ j^Q^ ^ ^ Q ^ Q-(^g^ spirals 

towards r = |y4| when xq — *■ — cxd and r = r+(s), ^ = ^"*'(s) reaches r = ri 
when Xq +00. Therefore r = \A\ is the boundary of a black hole. 
A similar results hold when B <Q. The difference is that the spiraling of 
the solutions towards r = 1^41 changes from the clockwise to the 

counter-clockwise or vice versa. 
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Example 4.2 Consider a generalization of Example 4.1 when 

v{x) = A{r)r + B{r)9, 

where ri < r < Tq, A{r),B{r) are smooth, 
A'^{ro) + B^{ro) = 1, ^^(r) + ^^(r) > 1 on [ri, tq), 5(r) > on 
[ri,ro], yl(r) + 1 has simple zeros ai, ...,arrti on (ri,ro), A(r) — 1 has simple 
zeros ■■■,Pm2 on ('^i; '^o), 7^ ttj;Vj, A;, |A(ri)| > 1. Equation (14.311) has 

the form 



(4.36) (^^(r) - 1)5,^ + {A{r)B{r) ± VA^ + B^ - 1)S^ = 0. 



Let 



(4.37) 



dr^(^)_A(\^ de+ _A{r)B{r) + VWTB^^ 
' A(r) + 1 ' 



(4.38) ^^A(.) + l, d^^Air)Bir) VWTW^_ 

^ ^ ds ^ ' ds A{r) - 1 

Note that the right hand sides in (I4.37p . (I4.38P are smooth since the 
singularities at A{r) + 1 = and A{r) — 1 = are removable. 

It follows from K37\\ . KM that 
{\x\ = aj} X R, J = 1, mi, = (3k} x R, A; = 1, m2, are boundaries 
of black or white holes. In particular if ai = mmj^k{c(j, Pk) and A{ri) < — 1 
then r = «! is the boundaries of a black hole. If Pi = minj^fc(«j, (3k) and 
A{ri) > 1 then r = /5i is the boundary of a white hole. 

Remark 4.4 (Axially symmetric metrics) 
Consider an equation (II. ip in x R where is a three-dimensional 
domain. Let {r,6,(p) be the spherical coordinates of x = (xi, 0:2, X3), i.e. 

Xi = r sin 6 cos ip, X2 = r sin 9 sin if, X3 = r cos 6. 
Suppose that g^'' are independent of ip : g^^ = g^^{r, 6). Consider 
characteristic surface S independent of xq and (/?: 



as _dS dr as 89 ^ _ ^ 2 3 
axk ar dxk do axk ' 



where 
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Substituting fl4.40p into fl4.39p we get an equation: 
(4.41) «"(,«)(f)%2a.^(,.)f § + (1)^0, 

We assume that a^^{r, 9) are also independent of (p. We consider (14.411) in 
two-dimensional domain u such that 5i < r < 52-, < ^3 < ^ < tt — ^4 when 
(r, 9) G uj. Imposing on u and a*-'(r, 9), 1 < ^, j < 2, the conditions of 
Theorems 13.11 and 14.11 we shall prove the existence of black or white holes in 
n X R with the boundary of the form 5*0 x S"^ x R where ip E , xq G R 

and 5*0 is a closed Jordan curve in the (r, 9) plane. 
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